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The  present  paper  considers  the  effect  of  stress  concentrations  on  the  damage  evolution  in  brittle  solids  with  disordered 
microstructure.  This  problem  is  herein  approached  studying  discrete  models  provided  with  a  local  load  sharing  rule.  Two  of 
these  models:  tight  bundle  parallel  bar  model  and  lattices,  are  employed  to  establish  universal  trends  in  brittle  cooperative 
deformation  processes. 


1.  Introduction 

The  stress  concentration  attributed  to  the 
presence  of  microcracks  and  microvoids  is  com¬ 
monly  regarded  as  one  of  the  leading  causes  of 
failure  in  materials.  As  a  result  of  the  defect 
induced  stress  concentration  material  adjacent  to 
the  existing  microdefect  is  much  more  likely  to 
fail  than  the  undamaged  regions  far  from  the 
defect  subjected  to  average,  or  far  field  stresses. 
Consequently,  a  rigorous  analysis  of  the  coopera¬ 
tive  brittle  phenomena  must  consider  spatial  dis¬ 
tributions  of  strengths  and  stresses.  In  contrast  to 
the  loose  bundle  parallel  bar  model  (considered 
in  the  first  part  of  this  study,  Krajcinovic  et  al. 
(1993))  it  is,  thus,  necessary  to  consider  not  only 
the  statistical  distribution  of  rupture  strengths 
but  also  the  preci.se  locations  of  weak  links  within 
the  system.  This  problem  is  obviously  much  more 
difficult  since  it  involves  a  more  complete  de¬ 
scription  of  the  microstructural  disorder. 
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The  proposed  strategy  is  to  extract  simple, 
universal  trends  and  features  of  the  macro  re¬ 
sponse  through  considerations  of  simplified  net¬ 
work  models  emphasizing  disordered  microstruc¬ 
tures.  The  universality  should  be  understood  as  a 
property  invariant  of  the  details  (higher  statistical 
momenta)  of  the  disorder,  damage  evolution,  etc. 
Once  determined,  universal  parameters  and 
trends,  being  common  to  an  entire  class  of  mate¬ 
rials  and  phenomena,  will  greatly  facilitate  formu¬ 
lations  of  more  sophisticated  models  necessary 
for  application  in  engineering  design. 


2.  Parallel  bar  model 

In  general,  a  local  load  sharing  rule  may  be 
introduced  rigorously  through  a  more  sophisti¬ 
cated  discretization  or  in  a  simpler,  but  more 
arbitrary,  manner  using  a  modified  version  of  the 
parallel  bar  model  allowing  for  unequal  distribu¬ 
tion  of  loads  released  by  the  ruptured  link.  To 
establish  connection  with  the  preceding  study  it 
seems  reasonable  to  consider  first  the  latter  alter¬ 
native.  In  a  parallel  bar  model  a  local  load  shar¬ 
ing  rule  may  be  postulated  a  priori  (directly)  or  by 
means  of  a  hierarchical  grouping  of  links  into 
cells. 
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2. 1.  Parallel  bar  model  with  local  load  sharing 

I'luctuations  of  forces  within  a  parallel  bar 
model  can  be,  for  example,  estimated  using  an  a 
priori  selected  load  sharing  rule.  For  example, 
Harlow  and  Phoenix  (1982)  assumed  that  the 
force  relea.sed  by  a  rupturing  link  is  equally  shared 
by  its  neare.st  neighbors.  Even  with  such  a  simple 
approximation  of  the  reality  rigorous  analyses 
turn  out  to  be  all  but  trivial.  In  fact  Harlow  and 
Phoenix  (1982)  were  only  able  to  derive  a  conser¬ 
vative.  lower  bound  on  the  system  strength  as¬ 
suming  that  the  failure  will  occur  when  the  two 
adjacent  links  in  the  bundle  rupture.  This  as¬ 
sumption  ma\  be,  indeed,  warranted  in  the  case 
of  ver\'  homogeneous  materials.  In  most  struc¬ 
tural  materials,  with  a  larger  scatter  in  rupture 
strengths,  macro  failure  is  typically  preceded  by 
the  formation  of  much  larger  clusters  of  ruptured 
links.  Unfortunately,  "the  complexity  (of  these 
statistical  models)  grows  enormously  with  increas¬ 
ing  bundle  size  and  simplification  techniques  are 
by  no  means  obvious"  (Phoenix,  1978). 

2.2.  Hierarchical  models 

In  many  cases  natural  systems  and  materials, 
such  as  tendons  (Cassidy  et  al.,  1990),  are  intrinsi¬ 
cally  hierarchical.  In  other  cases,  the  self-similar¬ 
ity  on  several  scales  is  the  product  of  manufactur¬ 
ing  processes  (solidification)  or  non-equilibrium 
growth  (vapor  deposition).  Finally,  the  scale  in¬ 
variance  may  also  be  the  consequence  of  the 
stress  directed  dilution  of  common  materials  with 
disordered  microstructure.  Indeed,  measure¬ 
ments  indicate  that  the  failure  surfaces  caused  by 
fracture  (Mandelbrot  et  al.,  1984),  faulting  and 
microcrack  growth  and  coalescence  are  character¬ 
ized  by  fractal  geometry.  Hence,  the  study  of 
hierarchical  structures  may  be  important  both  for 
their  own  sake  and  as  an  approximation  of  the 
local  load  sharing  rule. 

A  hierarchical  structure  studied  by  Smalley  et 
al.  (1985)  is  topologically  a  Cayley  tree  shown  in 
Fig.  1.  The  original  parallel  bar  model  (bottom 
level  in  Fig.  1)  containing  N  links  is  partitioned 
into  N /2  cells,  each  containing  two  links.  The 
rupture  strengths  of  individual  links  are  dis- 


Fig.  I.  A  parallel  bar  model  in  form  of  a  Cayley  tree  of 
functionality  two. 


tributed  according  to  an  appropriately  selected 
rupture  strength  distribution  function  p{f,).  The 
spatial  distribution  of  links  with  unequal  strength 
is  typically  varied  to  study  the  effect  of  different 
realizations  of  the  same  statistics. 

Initially,  al!  links  carr>'  identical  loads  f  =  F/N, 
where  F  is  the  externally  applied  tensile  load. 
When  a  link  ruptures  the  force  it  carried  is  re¬ 
leased  and  quasi-statically  transferred  to  the  re¬ 
maining  link  in  the  cell.  A  cell  fails  only  if  both 
constituent  links  rupture.  Consequently,  the  fail¬ 
ure  of  a  cell  of  order  (r)  represents  the  rupture 
of  a  corresponding  link  of  the  order  (r  -I-  I)  (see 
Fig.  1).  In  this  sense,  the  hierarchical,  Cayley 
tree,  model  introduces  a  simplified  load  sharing 
rule.  More  realistic  load  sharing  rules  may  be 
formulated  by  increasing  the  number  of  links  in  a 
cell  (functionality)  at  the  expense  of  the  simplicity 
and  tractability. 

The  transfer  of  forces  within  the  Cayley  tree 
was  quantified  by  Smalley  et  al.  (1985)  intro'luc- 
ing  conditional  probabilities  defining  the  chance 
that  a  cell  will  fail  when  the  force  from  a  rup¬ 
tured  link  is  transferred  within  the  same  cell  to 
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the  extent  link.  Omitting  the  details  of  the  deriva¬ 
tion,  available  in  Smalley  et  al.  (1985)  and  Herr¬ 
mann  and  Roux  (1990),  it  suffices  for  present 
purposes  to  list  the  final  recursive  expression 
relating  cumulative  (Weibull)  probabilities  of  fail¬ 
ure  Pif,)  in  cells  of  two  sequent  orders  (r)  and 
(r+  1), 

^(.+  0  =  2P‘^’[l  -  (1  (1) 

In  this  particular  case,  the  cumulative  proba¬ 
bility  function  P(f>  /,).  i.c.,  the  maximum  force 
which  can  be  entrusted  to  the  considered  hierar¬ 
chical  system,  can  be  solved  from  (1)  by  iteration. 
Instead,  Smalley  et  al.  (1985)  applied  the  renor¬ 
malization  group  method  according  to  which  the 
failure  corresponds  to  the  bifurcation  (or  unsta¬ 
ble  fixed)  point  in  the  (P*''*,  /’•'■*")  space.  For 
the  Weibull  rupture  strength  distribution  (Eq. 
(2.11)  in  Part  1)  (with  the  shape  parameter  «  =  2) 
the  bifurcation  point  corresponding  to  the  recur¬ 
sive  equation  (1)  is  P*  =  0.2063.  For  P  <  P*  the 
solution  tends  to  p  =  0  (stable  state),  while  for 
P>P*  the  solution  of  (1)  tends  to  p  =  1  (rup¬ 
ture). 

The  corresponding  value  of  the  maximum  force 
is  then  =  0.4807  ^  Ku  0.6798  Ku„  (where 
11+  and  are  the  Weibull  scale  parameter  and 
the  system  displacement  at  maximum  force  F  = 
Pn,).  In  comparison,  the  loose  bundle  (equal  load 
sharing)  estimate  for  the  same  distribution  of 
rupture  strengths  can  be  computed  from  Eq. 
(2.13)  in  Part  I  as  =  e“'/^  Ku^  =  0.6065  Ku^. 
The  difference  of  11%  between  these  two  esti¬ 
mates  is  obviously  a  direct  result  of  the  selected 
load  sharing  rule  (i.e,,  cell  size).  As  expected  with 
the  increasing  cell  size  the  influence  of  the  stress 
concentration  is  diminished.  According  to  the 
numerical  computations  for  the  Cayley  tree  hav¬ 
ing  cells  containing  four  links  (approximating  a 
two-dimensional  parallel  bar  model)  =  0.4327 
Kx+=  0.6119  Ku^  which  is  only  1.3%  in  excess  of 
the  loose  bundle  estimate.  This  result  provides  a 
first  but  nevertheless  convincing  evidence  that 
the  direct  interaction  of  defects  has  a  second- 
order  effect  on  the  system  response  during  the 
precritical  stage  of  the  deformation  process  and 
the  magnitude  of  the  maximum  force  P„,  itself. 


3.  Lattices 

On  the  next  level  of  sophistication  a  solid  may 
be  partitioned  into  irregular  polyhedra  (“finite 
elements”)  referred  to  as  Voronoi  tessellation 
(froth)  or  Wigner-Seitz  cells.  Cellular  represen¬ 
tation  of  solids,  which  predates  discretizations 
such  as  finite  elements,  are  natural  and  common 
in  modeling  stochastically  nonuniform  solids  with 
disordered  microstructure.  Dual  to  this  division 
into  cells  is  the  Delauney  simplicial  graph  (net¬ 
work,  see  Zallen  (1983)  and  Ostoja-Starzewski 
and  Wang  (1989))  shown  in  Fig.  2.  The  sides  of 
the  original  Voronrri  polyhedra  (“finite  ele¬ 
ments”)  are  orthogonal  to  the  bonds  (sides)  of 
the  Delauney  network  (lattice).  The  nodes  of  the 
simplicial  Delauney  network  are  equidistant  from 
the  point  of  intersection  of  two  graphs.  The  at¬ 
tribute  “simplicial”  refers  to  the  fact  that  the 
Delauney  division  consists  of  simplest  polyhedra 
(triangles  in  two  and  tetrahedra  in  three  dimen¬ 
sions)  or  simplexes.  It  can  be  shown  that  the 
Voronoi  polyhedra  and  Delauney  simplexes  can 
fill  the  entire  space. 

In  a  manner  analogous  to  one  suggested  by 
Englman  et  al.  (1984)  it  will  be  assumed  that  the 
load  carrying  capability  of  a  Voronoi  “finite  ele¬ 
ment”  is  exceeded  in  a  given  direction  when  the 
bond  of  the  dual  Delauney  graph  in  that  direc¬ 
tion  ruptures.  In  this  sense  the  considered  dis¬ 
cretization  is  a  generalization  of  the  “local  ap¬ 
proach  to  fracture”  concept  introduced  by 


Fig.  2.  Voronoi  polyhedra  and  its  dual  Delauney  network 
(bold  lines). 
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Lcmaitre  (1986)  and  successfully  applied  in 
finite-element  computations  thereafter. 

Introducing  Delauney  network  a  two-dimen¬ 
sional  continuum  is  approximated  by  an  articu¬ 
lated  network  of  links  joined  at  nodes.  The  nodes 
may  form  a  periodic  lattice  of  regular  geometry 
or  an  irregular,  topologically  disordered,  De¬ 
launey  network.  Structurally,  a  lattice  or  a  net¬ 
work  may  be  either  a  truss  (central-force  lattice) 
or  a  grid  (frame).  The  links  of  a  truss  are  sub¬ 
jected  only  to  a.\ial  forces.  The  elements  of  a  grid 
can  support  both  axial  forces  and  bending  mo¬ 
ments  preventing  relative  rotations  of  links  in  a 
node.  With  respect  to  the  connectivity  range  a 
node  may  be  connected  only  to  the  nearest  neigh¬ 
bors  (nodes  forming  the  first  coordination  group 
-  simple  trusses)  or  nodes  of  several  coordination 
groups  (complex  trusses).  The  present  study  will 
concentrate  only  on  some  of  the  many  possible 
alternatives  obtained  combining  features  men¬ 
tioned  above.  In  particular  only  the  planar  lat¬ 
tices  having  perfect,  periodic,  geometry  will  be 
considered  in  the  sequel. 

The  lattice  models  to  be  discussed  within  this 
study  should  not  be  confused  with  the  geometri¬ 
cally  similar  models  used  in  molecular  dynamics 
problems.  Present  study  considers  lattices  only  as 
elementary  discretizations  of  solids  allowing  for 
simple  analyses.  In  the  case  of  molecular  dynam¬ 
ics  models  nodes  actually  represent  the  molecules 
or  atoms.  Thus,  the  intermolecular  (bond)  forces 
must  be  derived  from  the  interaction  energy  typi¬ 
cally  approximated  by  the  “Lenard-Jones  6-12 
potential”. 

3.1.  Percolation  models 

In  view  of  the  substantial  literature  on  the 
percolation  theory  (see  Stauffer,  1985,  etc.)  a 
short  precis  of  basic  aspects  will  suffice  to  render 
the  following  discussion  at  least  partially  self-con¬ 
tained.  Only  the  bond  percolation  problem  will 
be  considered  and  p  will  denote  the  probability 
that  a  link  is  present  (and  q=\-p  that  the 
bond  is  absent  or  ruptured).  In  its  initial  state  of 
quenched  disorder  a  lattice  is  usually  weakened 
by  few  missing  links  randomly  distributed  over 
the  entire  lattice.  The  lattice  is  subsequently  di¬ 


luted  by  sequential  removal  of  links  which  can  be 
either  random  or  directed  by  applied  loads  (pos¬ 
tulating  an  appropriate  rupture  criterion).  In  the 
course  of  a  dilution  process  the  number  of  miss¬ 
ing  links  grows  along  with  the  probability  of 
emergence  of  larger,  irregularly  shaped  clusters 
of  missing  links.  Coincidentally,  the  stiffness  of 
the  lattice  decreases.  At  the  elastic  percohiiion 
threshold  p=p^^  the  stiffness  of  the  lattice  is 
reduced  to  zero.  At  this  point  the  initial  .state 
characterized  by  a  disordered  distribution  of  de¬ 
fects  abruptly  makes  a  transition  into  a  state 
dominated  by  a  single,  highly  localized  cluster  of 
missing  links.  At  the  percolation  threshold  a  tran¬ 
sition  takes  place  from  defects  characterized  by 
short-range  connectivity  to  a  defect  cluster  exem¬ 
plifying  long-range  connectivity  (localization).  The 
infinite  (spanning)  cluster  emerging  at  p  =  p^^. 
transects  the  infinite  lattice.  The  percolation 
thresholds  p^^  available  in  the  literature  pertain 
to  lattices  of  infinite  sizes.  Numerical  simulations 
are,  however,  for  obvious  reasons  performed  on 
large  but  finite  lattices.  The  effect  of  the  finite 
size  is,  subsequently,  eliminated  by  finite-size 
scaling. 

The  considered  dilution  process  is  character¬ 
ized  by  a  continuous  and  gradual  change  of  the 
microstructure  (number  of  extent  links)  which  at 
the  elastic  percolation  threshold  results  in  an 
abrupt,  qualitative  change  of  the  response  re¬ 
flected  in  the  singularity  of  the  system  compli¬ 
ance.  A  phenomenon  of  this  type  is  commonly 
referred  to  as  the  phase  transition.  As  in  the  case 
of  the  parallel  bar  model  the  Gibbs  energy  re¬ 
mains  continuous  (Krajcinovic  et  al.,  1993)  while 
its  second  derivative  with  respect  to  the  force 
becomes  singular  at  p=p„.  Consequently,  the 
considered  critical  phenomenon  (change  in  the 
connectivity  range)  is  a  second-order  phase  tran¬ 
sition.  Finally,  in  the  proximity  of  the  percolation 
threshold  the  lattice  stiffness  (or  less  rigorously 
the  elastic  modulus  on  the  macro-scale)  changes 
slowly  with  the  proximity  parameter  {p  -p^)  be¬ 
fore  vanishing  at  p  =  p^.  Hence,  the  lattice  stiff¬ 
ness  is  a  reasonable  candidate  for  the  order  pa¬ 
rameter. 

The  most  remarkable  aspect  of  the  percolation 
theory  is  that  the  .systems  belonging  to  the  same 
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universality  class  have  identical  asymptotic  behav¬ 
ior  in  the  vicinity  of  the  critical  state.  More 
specifically,  the  systems  (phenomena)  forming  a 
universality  class  have  identical  percolation 
thresholds  and  scaling  laws  defining  the  rate 
of  change  of  the  transport  properties  in  the  prox¬ 
imity  of  the  percolation  threshold.  While  the  per¬ 
colation  threshold  depends  only  on  tin.  dimen¬ 
sionality  d  and  the  lattice  type  (microstructure), 
the  scaling  law  for  the  stiffness  and  other  param¬ 
eters  exhibiting  singular  behavior  as  p-^p^  de¬ 
pends  only  on  the  dimensionality  d. 

In  its  conventional  form  the  percolation  model 
recognizes  only  two  alternatives;  a  link  is  either 
missing  (or  ruptured)  with  a  probability  q  =  \  —  p, 
or  is  present  with  a  probability  p.  All  links  are 
assumed  to  have  identical  cross  sectional  area  A 
and  elastic  modulus  E.  The  elastic  energy  of  the 
lattice  is  (modifying  slightly  the  expression  from 
Hansen,  1990) 

y 

+  EiZL-%g„e}.„  (2) 

ij 

where  tr,  is  the  displacement  of  the  node  t,  a 
unit  vector  located  at  node  i  and  directed  to¬ 
wards  the  node  j  and  the  change  of  angle 
between  the  links  ij  and  ik.  Also,  /  is  the  modu¬ 
lus  of  inertia  and  L  the  length  of  a  link.  Finally, 
L~'g^j-  L~'  if  the  link  is  present  and  zero  oth¬ 
erwise. 

3.I.J.  Central-force  lattices  (trusses) 

The  simplest  models  ensue  neglecting  the  link 
bending  stiffness  (EI/L^)  and  reducing  expres¬ 
sion  (2)  to  the  first  of  the  two  terms  on  the 
right-hand  side.  The  mean-field  theory  estimate 
for  the  central-force  lattice  stiffness  and  the  per¬ 
colation  threshold  for  a  random  dilution  can  be 
determined  in  a  very  simple  manner.  If  N  is  the 
number  of  nodes,  the  number  of  degrees  of  free¬ 
dom  is  then  (Nd)  where  d  is  the  lattice  dimen¬ 
sionality.  If  z  denotes  the  coordination  number 
(i.e.,  number  of  links  at  each  node)  then  the  total 
number  of  constraints  is  izNp/2).  During  a  ran¬ 
dom  dilution  process  a  truss  becomes  statically 


determinate  when  the  number  of  constraints 
equals  the  number  of  degrees  of  freedom,  i.e., 
when 

(3) 

Removal  of  a  single  additional  link  beyond  the 
percolation  limit  (3)  transforms  the  truss  into  a 
mechanism  unable  to  support  externally  applied 
loads.  In  (3)  the  subscript  (cen)  refers  to  the 
elastic  percolation  threshold  for  a  central-force 
lattice. 

During  the  initial  phase  of  the  process  the 
number  density  p  of  extant  links  is  well  above  the 
percolation  limit  p^^„.  During  this  phase  of  the 
deformation  process  the  missing  links  (defects), 
.scattered  over  the  lattice,  are  far  from  each  other. 
Forces  in  the  links  are  almost  equal  to  corre¬ 
sponding  forces  in  the  pristine  lattice  rendering 
the  self-consistent  theory  applicable.  Conse¬ 
quently,  the  lattice  stiffness  is  during  this  phase  a 
linear  function  of  p.  Thus,  the  mean-field  esti¬ 
mate  for  the  tangent  stiffness  of  a  randomly  di¬ 
luted  central-force  lattice  is 

^  P-Pccu 

K  1  ~  Peen 

where  K  is  the  stiffness  of  the  undamaged  (pris¬ 
tine)  lattice  for  which  p  =  1. 

Numerical  simulations  by  Feng  and  Sen  (1984), 
Sahimi  and  Goddard  (1986),  Beale  and  Srolovitz 
(1988)  and  others  on  diluted  triangular  two-di¬ 
mensional  central-force  lattices  lead  to  several 
important  conclusions.  In  all  cases  the  linear 
self-consistent  region  (4)  persisted  to  surprisingly 
small  tangent  stiffnesses  (/fy  =  0.05/f ).  The 
cross-over  region  connecting  the  mean-field  and 
percolation  regimes  was,  consequently,  very  short. 
Within  numerical  errors  the  percolation  thresh¬ 
old  was  found  to  be  =  0.65  which  for  z  =  6 
(triangular  lattice)  is  closely  fitted  by  the  mean- 
field  estimate  (3)  according  to  which  p^^  = 
2(2)/6  =  2/3. 

A  complete  analysis  of  a  large  central-force 
triangular  lattice  having  14,700  bonds  was  per- 
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formed  by  Beale  and  Srolovitz  (1^88)  who  esti¬ 
mated  the  percolation  threshold  (i.e.,  the  density 
of  existing  bonds  at  the  instant  of  vanishing  macro 
stiffness)  as  p  =  0.65  ±  0.005.  The  force- 
displacement  curve  of  the  considered  lattice,  sub¬ 
jected  to  tensile  tractions  at  two  ends,  was  found 
in  form  of  a  number  of  sharp  spikes  separated  by 
stretches  during  which  the  displacement  increases 
at  zero  force.  This,  physically  meaningless,  re¬ 
sponse  is  attributable  to  the  loss  of  internal  stabil¬ 
ity  of  the  truss  and  subsequent  rigid  body  pivoting 
of  the  overconstrained  segments  of  the  truss 
needed  to  align  the  extant  links  with  the  exter¬ 
nally  applied  forces.  This  instability  is  the  conse¬ 
quence  of  the  vectorial  nature  of  the  equilibrium 
equations  for  a  lattice  node  and  is,  in  fact,  not 
present  in  conduction  phenomena. 

Literature  focused  on  the  elastic  percolation 
problems  is  vague  on  a  rather  crucial  issue.  Al¬ 
most  as  a  rule  the  overall  (lattice  or  specimen) 
stiffness  is  confused  with  the  elastic  modulus. 
Even  more  importantly  it  is  not  obvious  whether 
the  scaling  law  pertains  to  the  tangent  or  secant 
stiffness  and  whether  the  percolation  threshold  is 
defined  as  a  dilution  level  at  which  the  tangent 
stiffness  Kj  or  the  secant  stiffness  K  reduces  to 
zero.  All  experiments  (Benguigui,  1984;  Sieradzki 
and  Li.  1986;  Benguigui  et  al.,  1987)  reported  in 
the  literature  w'ere  load  controlled.  Since  the 
post-peak  (softening)  segment  of  the  force-dis¬ 
placement  curve  is  unstable  a  load  controlled  test 
cannot  distinguish  between  K  =  0  and  Kj  =  0. 

The  arguments  in  the  first  part  of  this  study 
(Krajcinovic  et  al.,.  1993)  indicate  that  the  phase 
transition  occurs  al  the  apex  of  the  force-dis¬ 
placement  curve,  i.e.,  when  K-^  =  0.  However,  the 
infinite  cluster  emerges  at  this  point  only  in  force 
controlled  cases.  On  purely  physical  grounds  it  is 
appealing  to  conjecture  that  the  onsets  of  change 
in  connectivity  (percolation  threshold),  localiza¬ 
tion,  loss  of  homogeneity  in  displacement  field 
and  loss  of  ellipticity  of  the  governing  equations 
are  different  manifestations  of  the  same  phe¬ 
nomenon  (criticality).  However,  this  issue  can  be 
settled  conclusively  only  by  precise  force  and 
displacement  controlled  tests  which,  to  the  best 
knowledge  of  the  authors,  were  not  as  yet  re¬ 
ported  in  the  literature. 


3.L2.  Frames  (gruls) 

Two  major  problems  associated  with  the  cen¬ 
tral-force  lattice  model  are:  (a)  occurrence  of  the 
model  incured  instabilities;  and  (b)  the  fact  that 
the  elastic  percolation  threshold  significantly 
exceeds  the  percolation  threshold  for  the  con¬ 
ductivity  problems.  The  physics  and  mathematical 
structure  of  the  conductivity  and  elastic  percola¬ 
tion  problems  are  by  no  means  identical.  The 
Kirchhoff  equations  defining  the  partition  of 
electric  current  in  a  node  are  scalar  in  contrast  to 
the  nodal  equilibrium  equations  which  are  vecto¬ 
rial.  However,  in  two-dimensional  problems  two 
percolation  thresholds  should  coincide  since  the 
cluster  of  missing  links  traversing  the  lattice  (/\. ) 
transects  it  as  well  (p^^.)  (Krajcinovic  and  Basista, 
1991).  To  illustrate  the  role  of  the  internal  insta¬ 
bilities  it  suffices  to  examine  a  square  or  cubic 
lattice.  Since  neither  of  these  two  lattices,  even 
when  all  links  are  present,  has  any  shear  stiffness 
their  percolation  threshold  is  =  1  (as  ob¬ 
served  in  Roux  and  Guyon,  1985). 

There  are  at  least  two  ways  in  which  the  situa¬ 
tion  can  be  remedied.  Firstly,  the  local  instabili¬ 
ties  can  be  prevented  constraining  relative  rota¬ 
tions  of  links  in  nodes  assuming  built-in  condi¬ 
tions  at  each  node  and  endowing  the  links  with 
bending  rigidity.  Alternatively,  each  node  can  be 
connected  to  nodes  belonging  to  more  than  one 
coordination  group  to  reduce  the  chance  for  the 
occurrence  of  internal  instabilities.  At  this  point 
only  the  first  of  these  two  alternatives  will  be 
examined. 

Determination  of  the  elastic  peicolation 
threshold  p^^  for  a  grid  necessitates  extensive 
numerical  simulations.  According  to  the  results  of 
simulations  reported  in  Kantor  and  Webman 
(1984),  Feng  et  al.  (1985).  Sahimi  (1986)  and 
others,  it  turns  out  that  for  a  grid  p^^  =  p^.  Thus, 
as  expected  on  purely  geometric  grounds,  elastic 
and  conductivity  percolation  thresholds  (once  the 
internal  instabilities  are  eliminated)  are  indeed 
identical.  This  facilitates  application  of  the  perco¬ 
lation  theory  in  mechanics  since  the  conductivity 
percolation  thresholds  for  several  different  peri¬ 
odic  lattices  are  available  in  Stauffer  (1985)  and 
other  standard  monograph  on  the  percolation 
theory.  The  effect  of  the  rotational  constraints  in 
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lattice  nodes  is  very  significant.  For  example,  the 
fraction  of  bonds  needed  to  maintain  the  load 
carrying  capability  is  lowered  from  0.645  to  0.347 
for  triangular  lattice  and  from  1.0  to  0.5  for  the 
quadratic  lattice. 

Finally,  an  argument  can  be  made  that  the 
regularity  of  the  considered  lattices  has  no  effect 
on  the  conclusions.  Being  simplicial,  a  planar 
Delauney  network  has  an  identical  coordination 
number  z  =  6  as  the  regular  triangular  lattices. 
Since  the  product  zp^  ~  2.0  ±  0.2  (where  the  su¬ 
perscript  b  stands  for  bond  percolation)  is  a  di¬ 
mensional  invariant  for  all  examined  lattices  (Zal- 
len,  1983)  it  seems  reasonable  to  conclude  that 
the  bond  percolation  thresholds  for  triangular 
and  Delauney  lattices  are  similar  in  magnitude. 
Consequently,  the  pre-peak  (hardening)  response 
of  a  triangular  lattice  is  not  affected  by  the  regu¬ 
larity  (periodicity)  of  the  lattice.  An  analogous 
conclusion  was,  in  the  context  of  the  sol-gel 
transition  reached  by  de  Gennes  (1979,  Ch.  V.2.8). 

3.1.3.  Directed  percolation 

The  above  listed  data  are  related  to  perfectly 
random  dilution  of  lattices.  In  the  course  of  ran¬ 
dom  dilution  processes  the  lattice  stiffness  tensor 
changes  isotropically.  A  glance  at  the  lattices 
studied  by  Beale  and  Srolovitz  (1988)  indicates 
that  in  the  case  of  dilution  processes  driven  by 
arbitrarily  oriented  externally  applied  forces  (un¬ 
iaxial  in  their  case)  only  some  of  the  components 
of  the  stiffness  tensors  become  singular  at  the 
percolation  threshold.  The  links  aligned  with  the 
externally  applied  tractions  are  subjected  to  the 
largest  forces  and  are,  therefore,  much  more  likely 
to  fail  than  other  links. 

In  the  case  of  directed  percolation  infinite 
cluster  of  missing  links  forms  in  a  direction 
roughly  perpendicular  to  the  applied  tensile  trac¬ 
tions.  Therefore,  it  requires  rupture  of  fewer 
links  than  necessary  for  random  (isotropic)  dilu¬ 
tion.  According  to  Kinzel  (1983)  the  directed 
bond  percolation  thresholds  are:  p^^  =  0.479  for  a 
triangular  lattice  and  0.6445  for  a  square  lattice. 
Therefore,  the  Beale  and  Srolovitz  (1988)  esti¬ 
mate  of  =  0.65  should  be  compared  with 
Peed  ==  0.479  for  directed  percolation  and  not,  as 
typically  assumed,  to  =  0.347  for  the  isotropic 


(random  dilution)  case.  The  difference  between 
/>„„  and  is  attributable  to  the  direction  of 
applied  forces,  while  the  rest  (to  p^^)  is  due  to 
the  lack  of  rotational  constraints  in  nodes  (i.e., 
internal,  model  induced,  instabilities). 


4.  Distributed  strength  lattices 

One  of  the  basic  shortcomings  of  conventional 
percolation  models  is  inherent  to  the  bimodal 
distribution  of  the  rupture  strengths  of  links  which 
can  be  either  equal  to  zero  (initial  damage)  or  to 
some  positive  constant  /,,.  In  view  of  the  discus¬ 
sion  of  the  parallel  bar  models  in  the  first  part  of 
this  study  (Krajcinovic  et  al.,  1993)  it  is  obvious 
that  the  band-width  of  the  rupture  strength  distri¬ 
bution  has  a  strong  influence  on  the  type  of  the 
response  in  both  qualitative  and  quantitative 
sense.  However,  analyses  of  disordered  lattices 
assembled  from  links  having  arbitrary  distributed 
rupture  strengths  and  stiffness  require  large-scale 
numerical  simulations.  First  statistical  studies  of 
triangular  lattices  conducted  by  Sahimi  and  God¬ 
dard  (1986)  were  directed  towards  investigation 
of  the  influence  of  different  link  stiffnesses  and 
different  link  rupture  thresholds  on  the  lattice 
response.  These  analyses  demonstrated  that  by 
changing  the  distribution  of  link  rupture  strengths 
it  becomes  possible  to  change  the  failure  mode 
from  one  in  which  the  spanning  crack  is  of  very 
regular  geometry  to  one  characterized  by  a  span¬ 
ning  crack  of  very  complex  (fractal)  geometry.  In 
the  first  case,  characterized  by  a  narrow  link 
strength  bandwidth,  the  failure  mode  is  brittle 
since  almost  all  of  the  energy  is  dissipated  to 
form  the  spanning  crack.  The  failure  in  this  rela¬ 
tively  homogeneous  material  can  be  attributed 
solely  to  the  stress  concentration  of  a  single  de¬ 
fect  of  preferential  geometry.  In  the  latter  case, 
characterized  by  a  more  substantial  link  strength 
bandwidth,  a  large  number  of  links  is  ruptured 
without  being  the  part  of  the  spanning  cluster 
and  promoting  the  macro-failure.  In  other  words, 
in  this  damage  tolerant  heterogeneous  solid  fail¬ 
ure  is  preceded  by  significant  dissipation  of  en¬ 
ergy  and  attendant  accumulation  of  damage. 

A  truly  comprehensive  study  of  triangular  cen- 
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tra!-force  lattices  consisting  of  links  with  identical 
stiffnesses  and  randomly  distributed  rupture 
strengths  was  recently  completed  by  Hansen  ct  al. 
(1989).  At  each  end  the  lattice  was  supplied  by  a 
rigid  member  (bus)  ensuring  identical  displace¬ 
ment  of  all  nodes  it  connects.  Periodicity  condi¬ 
tions  were  enforced  in  lateral  direction  to  elimi¬ 
nate  the  end  effects.  The  quenched  disorder  was 
introducing  assuming  uniform  distribution  of  link 
rupture  strengths  pif,)  =  constant  (modeling  the 
initial  damage  by  links  of  zero  strength).  On  such 
a  lattice  Hansen  et  al.  ( 1989)  performed  repeated 
computations  for  a  large  number  of  different 
physical  realizations  of  the  same  statistics.  The 
only  difference  between  individual  simulations 
consisted  of  different  spatial  distribution  of  links 
strengths.  To  ascertain  the  effect  of  the  lattice 
size  Hansen  et  al.  (1989)  considered  four  differ¬ 
ent  lattice  sizes  (with  A  =  4.  8,  16  and  32  number 
of  rows,  see  Fig.  3).  Computations  for  the  lattice 
of  smallest  size  (A  =4)  were  considered  statisti¬ 
cally  unreliable  since  the  individual  defects  (mis¬ 
sing  links)  occupy  a  large  portion  of  the  lattice 
from  the  very  beginning  of  the  dilution  process. 

Computations  were  performed  following  the 
rules  of  conventional  truss  analy.ses.  Forces  in 
each  link  were  computed  for  each  external  dis¬ 
placement  increment.  A  link  was  removed  and  its 
force  redistributed  to  extant  links,  in  accordance 
with  the  equilibrium  conditions,  whenever  the 
force  in  this  link  exceeded  its  strength  f,.  Hansen 
et  al.  (1989)  considered  a  symmetric  rupture  con¬ 
dition  making  no  difference  between  tension  and 
compicssion.  However,  the  ensuing  departure 
from  physical  reality  of  fracture  was  negligible 


since  the  considered  tensile  loading  was  not  con¬ 
ducive  to  emergence  of  large  compressive  forces. 
Hence,  since  the  lattice  geometry'  was  periodic, 
and  analysis  deterministic,  the  disorder  was  a 
direct  consequence  of  the  nondeterministic  distri¬ 
bution  of  the  initial  damage  and  rupture  strengths 
(quenched  disorder). 

The  macro-variables  defining  the  kinematics  of 
the  considered  dissipative  process  are  the  dis¬ 
placement  u  of  the  lattice  (change  of  the  distance 
separating  the  rigid  members)  and  an  appropri¬ 
ately  defined  damage  variable  D.  On  micro-scale 
the  history'  is  recorded  by  the  number  of  broken 
links  n.  Since  n.  or  D.  are  the  only  history 
recording  variables  all  macro-variables  and  trans¬ 
port  properties  (such  as  lattice  stiffness)  mu.st  be 
defined  as  a  function  of  D  and  all  micro-variables 
as  a  function  of  n.  To  determine  the  statistics  of 
the  proce.ss  the  results  of  computations  arc  aver¬ 
aged  over  the  entire  ensemble  of  physical  realiza¬ 
tions  selecting  the  number  n  of  ruptured  links  as 
the  control  variable.  In  view  of  the  finite  size  of 
the  lattices  the  force-displacement  relation  w'as 
sought  by  Hansen  et  al.  (1989)  in  form  of  the 
following  finite-scaling  Ansatz 

F  =  K{\-^l3au)u.  (5) 

The  lattice  stiffness  K  was,  for  convenience,  taken 
to  be  equal  to  unity  by  Hansen  et  al.  (1989).  The 
parameter  a  was  then  selected  to  ensure  the  best 
fit  between  the  Ansatz  (5)  and  the  computational 
results  over  the  entire  spectrum  of  parameters 
(lattice  sizes).  It  was  not  entirely  surprising  that 
the  parameter  a  itself  was  found  to  be  a  function 
of  lattice  size  L.  On  the  basis  of  these  computa¬ 
tions  Krajcinovic  and  Basista  (1991)  proposed  a 
simple  equality 

a=/3A'-^,  (6) 

which  for  p  =  3/4  (as  suggested  in  Hansen  et  al., 
1989)  fits  the  computed  data  exceedingly  well. 

Having  the  above  expression  (6)  for  the  param¬ 
eter  a  the  force-displacement  relation  can  be 
cast  in  the  familiar  form  (see  part  1,  Eq.  (2.2)) 
(allowing  additionally  for  the  lattice  size  effect) 

F(A)=/w[l-ra(A)]«.  (7) 
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The  damage  parameter  a>  can  be  shown  (for 
the  selected  uniform  distribution  of  link  rupture 
strengths)  to  be  linearly  proportional  to  the  dis¬ 
placement  u 


The  damage  parameter  w  can  be  also  rewritten 
as  a  function  of  the  externally  applied  force  F  in 
form  of 

a)(A,  F)  =U.5{l  -  [l  -4v^/3A-''{F/A:  )]'  ■}. 


where  as  before  /t)  =  3 /4. 

Similarly,  the  displacement  corresponding  to 
the  apex  of  the  force-displacement  curve  is  from 
(7)  and  (8) 

=  (10) 

where  the  subscript  m  denotes  the  value  of  the 
variable  corresponding  to  the  maximum  force 
to  which  the  lattice  can  be  subjected  (rupture 
force  in  a  load  controlled  test).  The  number  of 
ruptured  links  at  the  apex  of  the  force-displace¬ 
ment  curve  is 

n„  =  A^*'/2)3.  (11) 

Substitution  of  (10)  into  (8)  and  (7)  leads  to 
analytical  expressions  for  the  damage  at  the  apex 
and  maximum  force  which  are  both  independent 
of  the  lattice  size 

a.„  =  0.5  and  F„  =  0.5  Ku^.  (12) 

From  (7)  the  secant  modulus  K  =  Kil  -  m)  at  the 
apex  is  according  to  (12)  equal  to  0.5  K.  The 
analytical  expression  for  the  number  of  ruptured 
bonds  at  the  apex  is 

n,„  =  u„\'j5.  (13) 


The  self-consistent  estimate  for  the  displace¬ 
ment  can  be  readily  derived  in  a  manner  indi¬ 
cated  in  Hansen  et  al.  (1089).  For  the  considered 
tensile  loading  only  the  diagonal  links  are  likely 
to  rupture.  Since  the  total  number  of  diagonal 
links  is  2A',  the  damage  D  (number  density  of 
ruptured  links)  is.  according  to  Eq.  (2.3)  of  Part  1. 
equal  to  n/lX'.  Thus,  the  average  strain  is 

<e>  =  e^^.  (14) 

where  'S  the  strain  at  which  the  strongest  link 
ruplure.s.  i.e..  for  which  D  =  \. 

At  the  same  time  since  iii/X)  is  the  displace¬ 
ment  of  a  lattice  node 


where  the  length  of  the  diagonal  link  /  is  equal  to 
unity.  Hence,  the  expression  (13)  may  indeed  be 
derived  from  (14)  and  (15)  as  a  self-consistent 
estimate  of  the  number  of  ruptured  links  at  the 
state  corresponding  to  the  apex  of  the  force -dis¬ 
placement  curve. 

The  parameter  a,  number  and  density  of  rup¬ 
tured  links  displacement  force  F„  and 
damage  at  the  apex  of  the  force-displace¬ 
ment  curve  (failure  in  the  stress  controlled  test) 
were  computed  from  the  above  expressions  and 
arranged  in  Table  1.  The  corresponding  numeri¬ 
cal  data  read  off  Fig.  6  in  Hansen  et  al.  (1989)  are 
added  in  parentheses.  As  already  mentioned  the 
numerical  results  for  A  =  4  are  omitted  as  sug¬ 
gested  in  Hansen  et  al.  (1989). 

The  remarkable  accuracy  with  which  the  above 
expressions,  derived  using  the  simple  relationship 
(6)  for  the  parameter  a,  fit  the  numerical  simula¬ 
tions  provides  a  strong  indication  that  the  lattice 


Table  1 

Parameter  a  (Eq.  (5)),  number  and  density  of  ruptured  bonds,  displacement,  force,  and  damage  at  the  apex  of  the  force-displace¬ 
ment  curve 
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response  during  the  liardening  phase  of  the  de- 
formatioii  proeess  (i.e.,  along  the  ascending  part 
of  the  force  displacement  cur\e)  depends  en¬ 
tirely  on  the  volume  averages  of  the  disorder 
(spatial  distribution  of  ruptured  links).  This  is, 
indeed,  the  reason  behind  the  close  agreement 
between  the  self-consistent  estimates  and  the  re¬ 
sults  of  numerical  computations  from  Hansen  et 
al.  (198*)).  Moreover,  since  the  analytical  esti¬ 
mates  of  the  pre-peak  (hardening)  .segment  of  the 
response  obtained  by  the  parallel  bar  model  are 
in  close  agreement  with  the  numerical  simula¬ 
tions  performed  on  lattices  the  latter  models 
should  be  used  only  when  the  emphasis  is  on  the 
post-peak  (.softening)  behavior. 

According  to  computations  in  Hansen  et  al. 
(1989)  both  the  total  number  of  ruptured  bonds 
/ii„  and  the  density  of  the  ruptured  bonds  (/i,,,//^ ) 
at  the  apex  strongly  depend  on  the  specimen  size. 
However,  the  magnitude  of  the  damage  variable 
at  the  apex  depends  only  on  the  selected  dist  ibu- 
tion  of  ruptured  strengths  as  already  ascertained 
by  the  expressions  derived  in  the  preceding  part 
of  this  paper  (Krajcinovic  et  al.,  1993)  for  a  paral¬ 
lel  bar  system.  It  is  important  to  notice  that  in 
contrast  to  the  parallel  bar  model  it  becomes 
necessary  to  make  distinction  between  the  density 
of  ruptured  bonds  D  =  (n/N)  and  the  parameter 
ti>  defining  the  reduction  of  the  secant  modulus 
K  =  K{\-u)).  Both  of  these  parameters,  as  ex¬ 
pected,  depend  on  the  microstructure  (selected 
distribution  of  the  link  rupture  strengths  p(/,)). 
However,  only  the  parameter  oj,  or  more  accu¬ 
rately  the  magnitude  of  the  secant  modulus  A'(ai), 
is  invariant  of  the  lattice  size.  Since  the  magni¬ 
tude  of  the  specimen  (lattice)  secant  modulus 
Kiu))  is  readily  measurable  in  experiments  this 
conclusion  is  potentially  very  important  for  the 
selection  of  the  internal  variable  quantifying  the 
damage. 

The  other  important  .set  of  conclusions  in 
Hansen  et  al.  (1989)  concerns  the  post-peak  or 
softening  part  of  the  lattice  response.  As  indi¬ 
cated  above  the  ascending  seg  ents  of  the 
rescaled  force-displacement  relation  (7).  in  con¬ 
junction  with  (9).  determined  for  different  lattice 
sizes  A  were  found  to  collapse  on  a  single  master 
curve.  A  single  parameter  /3,  in  conjunction  with 


the  finite-size  scaling  lasv  (5)  and  (6).  suffices  to 
eliminate  the  inlluence  of  the  lattice  size  A,  de¬ 
tails  (higher  statistical  momenta)  of  the  damage 
distribution  and  microstructurc  (precise  location 
of  links  having  given  rupture  strengths  within  the 
lattice).  The  scatter  of  results  about  their  means 
was  found  to  be  very  small  as  well.  In  sui.imary'. 
the  response  of  the  lattice  along  the  ascending 
segment  of  the  force -displacement  curve,  includ¬ 
ing  the  apex  itself,  can  be  with  acceptable  accu¬ 
racy  estimated  using  the  scif-consistent  models 
(or  some  other  first-order  effective  continuum 
model). 

In  the  pre-peak,  effective  continuum  regime 
the  defects  (clusters  of  missing  links)  arc  small 
and  far  away  from  each  other.  The  slre.ss  fluctua¬ 
tions  arc  local  and  spread  over  a  small  part  of  the 
lattice.  Consequently,  the  sequence  of  the  link 
ruptures  is  dominated  by  the  distribution  ot  link 
strengths  (de  Arcangelis,  1990).  Within  this  phase 
of  the  deformation  process  spatial  distribution  of 
links  of  different  strength  has  a  second-order 
effect  on  the  response.  The  lattice  response  is 
local  and  governed  by  the  averages  of  the  in¬ 
volved  fields. 

The  modeling  of  the  .softening  response  turned 
out  to  be  everything  but  simple.  The  numerical 
results  of  Hansen  et  al.  (1989)  strongly  suggest 
development  of  a  multi-fractal  force  distribution 
at  the  incipient  rupture  of  the  lattice.  Within  the 
post-peak  regime,  characterized  by  significant  de¬ 
fect  densities,  the  stress  fluctuations  surrounding 
the  defects  spread  over  most  of  the  lattice. 
Therefore,  each  link  is  subjected  to  a  different 
force.  The  local  (rather  than  the  average)  stress 
becomes  a  dominant  factor  of  the  lattice  response 
near  the  point  of  rupture.  Omitting  the  details, 
available  in  Hansen  et  al.  (1989),  Herrmann  et  al. 
(1989),  de  Arcangelis  (1990)  and  Hansen  (1990), 
it  suffices  to  state  that  the  scaling  behavior  of  the 
lattice  at  the  incipient  rupture  requires  an  infinity 
of  exponents  defining  different  statistical  mo¬ 
menta  of  the  distribution  of  forces  in  cutting 
links.  The  post-peak  response  preceding  the  lat¬ 
tice  rupture  depends  on  the  higher  statistical 
momenta  of  the  initial  (quenched)  disorder,  i.e., 
unlikely  event.s.  The  scatter  of  the  results  in  the 
softening  regime  is  much  more  pronounced  and 
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quite  persistent.  Experimental  evidence  of  this 
phenomenon  is  well  documented  in  the  case  of 
heterogeneous  materials  such  as  concrete  (see, 
for  example,  Hegemier  and  Reed,  1985). 

Finally,  it  should  not  be  construed  that  the 
multi-fractal  nature  of  the  force  distribution  in 
links  during  the  softening  regime  has  anything  to 
do  with  the  discretization  or  modeling  of  the  solid 
by  a  truss.  Identical  analyses  performed  for  a  grid 
(frame)  by  Herrmann  et  al.  (1989)  demonstrated 
iden'ical  fends  leading  to  the  same  conclusions. 

5.  Summary  and  conclusions 

The  study  presented  in  this  two-part  paper  has 
a  qualitative  character.  The  considered  dis¬ 
cretizations  are  insufficiently  sophisticated  to 
provide  for  a  rigorous  and  quantitative  descrip¬ 
tion  of  the  behavior  of  brittle  solids  with  disor¬ 
dered  microstructure.  Additionally,  all  considera¬ 
tions  in  *' 's  study  were  based  on  analyses  of 
two-dimen.'ional  systems  overestimating  the  role 
of  the  defect  interactions  on  the  macro  response. 
Finally,  all  lattices  considered  in  this  study  were 
limited  in  size  emphasizing  the  size  effect.  Never¬ 
theless,  unencumbered  by  complex  mathematical 
structure  and  amenable  to  analytical  solutions 
and  inexpensive  numerical  simulations  the  con¬ 
sidered  models  arc  valuable  in  discerning  the 
dominant  and  universal  trends  in  the  behavior  of 
disordered  solids. 

The  macro  response  of  a  solid  with  disordered 
structure  during  a  brittle  deformation  process 
consists  of  three  distinctly  different  regimes: 

(a)  During  the  initial  regime,  characterized  by 
a  low  density  of  small  defects  scattered  through¬ 
out  the  solid,  the  macro  response  is  local.  Trans¬ 
port  parameters,  such  as  lattice  stiffness,  can  be 
computed  from  volume  averages  of  defect:  and 
the  size  effects  are  readily  determined  using  ap- 
prop  iate  finite-size  scaling.  Stress  fluctuations 
are  both  small  and  local.  Owing  to  its  effective 
continuum  nature  the  hardening  segments  of  the 
response  of  systems  containing  similar  densities 
of  defects  are  similar  as  well. 

(b)  The  character  of  the  response  changes 
drastically  in  the  vicinity  of  the  elastic  percolation 


threshold.  The  scaling  law  for  the  system  tangent 
stiffness  is  characterized  by  a  strongly  nonlinear 
dependence  on  the  proximity  parameter.  Yet,  the 
self-consistent  estimates  of  the  damage  and  se¬ 
cant  moduli  at  the  percolation  threshold  are  in 
close  agreement  with  numerical  simulations. 

(c)  The  response  in  the  post-peak  (softening) 
regime  depends  on  the  higher  (extreme)  statisti¬ 
cal  momenta  of  the  initial  disorder  and  stress 
distribution  leading  to  a  large  test-to-test  scatter. 
Consequently,  the  predictions  of  the  expected 
values  of  the  ultimate  rupture  strength  of  the 
system  are  of  limited  utility. 

From  the  modeling  viewpoint  it  is  important  to 
notice  that  the  effective  continuum  models  (Horii 
and  Nemat-Nasser,  1983;  Krajeinovic  and 
Sumarac,  1989;  Krajeinovic,  1989;  Nemat-Na.sser 
and  Hori,  1990.  etc.)  can  be  used  along  the  entire 
length  of  the  hardening  regime  including  the  predic¬ 
tion  of  the  percolation  threshold  (onset  of  localiza¬ 
tion).  If  further  substantiated  this  conclusion  may 
render  development  of  higher-order  effective 
continua  models  unnecessary.  Identifying  the 
considered  process  of  damage  evolution  as  a  sec¬ 
ond-order  phase  transition  justifies  application  of 
several  novel  methods  of  statistical  physics.  More 
importantly  it  becomes  possible  to  identify  uni¬ 
versal  parameters  and  trends  in  deformation. 
These  universal  parameters  are  invariant  w'ith 
respect  to  discretization  (lattice  regularity)  and 
details  of  the  microstructural  disorder.  They  are 
common  to  a  wide  spectrum  of  materials  and 
phenomena. 

At  this  point  it  seems  certain  that  the  secant 
stiffness  modulus  represents  a  proper  choice  for 
the  internal  variable  quantifying  accumulated 
damage.  In  the  pre-peak  response  the  ^ecant 
moduli  are  readily  measured  in  tests  (Lemaitre 
and  Chaboche,  1978)  in  addition  to  having  at¬ 
tributes  of  universality  (invariance  with  respect  to 
the  specimen!  size,  dilution  sequence,  details  of 
the  microstructural  disorder,  etc,).  In  three-di¬ 
mensional  problems  the  stiffness  tensor  inher¬ 
ently  incorporates  anisotropy  which  is  one  of  the 
criteria  defining  the  universality  group  (Sengers, 
1985).  Finally,  the  damage  effective  continua  mi¬ 
cromechanical  models  (Horii  and  Nemat-Nasser, 
1983;  Krajeinovic,  1989,  etc.)  is  also  measured  as 
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the  change  in  the  stiffness  tensor.  Hence,  the  two 
theories,  valid  for  different  ranges  of  microdefect 
densities,  are  mutually  compatible. 

In  conclusion,  the  present  two-part  study  of  a 
set  of  rather  primitive  discrete  models  leads  to 
some  interesting  conclusions  related  to  coopera¬ 
tive  phenomena  in  brittle  deformation.  The  study 
helps  to  classify  the  mechanical  damage  evolution 
models  as  scconcl-onlcr  phase  transitions  justifying 
application  of  novel  statistical  physics  theories. 
The  study  also  identifies  necessary  thermody¬ 
namic  variables  defining  the  state  and  the  change 
of  state  in  a  deformation  process  dominated  by 
microcrack  evolution.  Most  importantly,  the  pre¬ 
sent  study  establishes  unirersal  trends  and  param¬ 
eters  defining  critical  state  and  behavior  for  the 
considered  class  of  problems. 
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